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Abstract 

We  consider  kinematic  chains  evolving  on  a  finite-dimensional  Lie  group  G  under 
nonholonomic  constraints,  where  snake-like  global  motion  is  induced  by  shape  variations 
of  the  system.  In  particular,  we  consider  the  case  when  the  evolution  of  the  system  is 
restricted  to  a  subspace  h  of  the  corresponding  Lie  algebra  G,  where  h  is  not  a  subalgebra 
of  Q  and  it  can  generate  the  whole  algebra  under  Lie  bracketing.  Such  systems  are 
referred  to  as  G-snakes.  Away  from  certain  singular  configurations  of  the  system,  the 
constraints  specify  a  (partial)  connection  on  a  principal  fiber  bundle,  which  in  turn  gives 
rise  to  a  geometric  phase  under  periodic  shape  variations.  This  geometric  structure 
can  be  exploited  in  order  to  solve  the  nonholonomic  motion  planning  problem  for  such 
systems. 

G-snakes  generalize  the  concept  of  nonholonomic  Variable  Geometry  Truss  assem¬ 
blies,  which  are  kinematic  chains  evolving  on  the  Special  Euclidean  group  SE( 2)  under 
nonholonomic  constraints  imposed  by  idler  wheels.  We  examine  in  detail  the  cases  of 
3-dimensional  groups  with  real  non- abelian  Lie  algebras  such  as  the  Heisenberg  group 
H( 3),  the  Special  Orthogonal  group  50(3)  and  the  Special  Linear  group  SL( 2). 
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1  Introduction 

Of  significant  interest  among  mechanical  systems  subject  to  nonholonomic  constraints  are 
those  wherein  variations  of  shape  induce,  under  the  influence  of  the  constraints,  a  global  motion 
of  the  system.  A  well-known  example  is  that  of  a  free-floating  multibody  system  in  space 
(e.g.  robotic  manipulators  mounted  on  orbiting  satellites),  where  periodic  movements  of  the 
joints  induce  a  reorientation  of  the  system,  under  the  nonholonomic  constraint  of  conservation 
of  angular  momentum  (Krishnaprasad  [1990];  Marsden,  Montgomery  &  Ratiu  [1990]). 

Inspired  by  the  experimental  work  of  Joel  Burdick  and  his  students  at  Caltech  (Chirikjian  & 
Burdick  [1991];  Chirikjian  &  Burdick  [1993]),  a  novel  system  that  uses  the  above  principle  for  land 
locomotion  was  introduced  in  (Krishnaprasad  &  Tsakiris  [1994]).  There,  a  Variable  Geometry 
Truss  (VGT)  assembly  consisting  of  longitudinal  repetition  of  truss  modules,  each  one  of  which 
is  equipped  with  idler  wheels  and  linear  actuators  in  a  planar  parallel  manipulator  configuration, 
uses  periodic  changes  of  the  shape  of  each  module  to  produce  global  motion.  The  locomotion 
principle  is  not  based  on  direct  actuation  of  wheels,  but  rather  on  the  nonholonomic  constraints 
imposed  on  the  motion  of  the  system  by  the  rolling  without  slipping  of  the  idler  wheels  of  each 
module  on  the  supporting  plane.  This  results  in  a  snake-like  motion  of  the  VGT  assembly, 
which  is  not  too  far,  at  least  in  principle,  from  certain  modes  of  actual  snake  locomotion  (Hirose 
[1993]).  Both  the  shape  and  the  configuration  of  the  VGT  assembly  can  be  described  by  elements 
of  the  Special  Euclidean  group  SE( 2),  the  group  of  rigid  motions  on  the  plane.  A  system  like 
the  VGT  assembly  constitutes  a  kinematic  chain  evolving  on  this  matrix  Lie  group,  with  the 
corresponding  velocities  given  by  elements  of  the  Lie  algebra  of  SE( 2).  Of  these  velocities,  the 
shape  variations  can  be  considered  as  the  controls  of  the  system  and  they  are  referred  to  as 
shape  controls.  The  nonholonomic  constraints  allow  us  to  express  the  global  motion  of  the  VGT 
assembly  as  a  function  of  the  shape  and  of  the  shape  controls  and  to  formulate  motion  control 
strategies  under  periodic  shape  controls. 

This  situation  can  be  generalized  to  kinematic  chains  evolving  on  an  arbitrary  (matrix)  Lie 
group  G  under  a  certain  class  of  nonholonomic  constraints.  In  particular,  we  are  interested  in 
groups  with  a  real  finite-dimensional  non-abelian  Lie  algebra  Q  (of  dimension  n)  and  {l—  1)- 
module  kinematic  chains  evolving  on  them,  where  the  constraints  force  the  velocities  of  the  system 
to  lie  in  a  subspace  of  G,  which  is  not  a  subalgebra  of  G  but  which  generates  the  whole  algebra  G 
under  Lie  bracketing.  We  refer  to  systems  of  this  type  as  G-snakes  and  observe  that  they  possess 
an  interesting  geometric  structure:  When  l  =  n  and  the  codimension  of  the  constraints  is  one, 
the  configuration  and  shape  spaces  of  the  system  specify  a  principal  fiber  bundle  (Bleecker  [1981]; 
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Nomizu  [1956])  and  the  nonholonomic  constraints  specify  a  (partial)  connection  on  it,  at  least 
away  from  certain  configurations  which  we  call  nonholonomic  singularities  (higher  codimension 
cases  will  be  treated  elsewhere). 

In  section  2  of  this  paper,  we  consider  an  ( l  —  l)-module  kinematic  chain  evolving  on  an 
n-dimensional  Lie  group.  The  Wei-Norman  representation  of  G  (Wei  &  Norman  [1964]),  which 
expresses  each  element  of  the  group  as  a  product  of  the  one-parameter  subgroups  of  6',  and  the 
notion  of  the  adjoint  action  of  G  on  Q  allow  us  to  express  in  a  compact  form  how  the  motion  of 
each  module  of  the  kinematic  chain  relates  to  that  of  the  other  modules  and  to  the  global  motion 
of  the  system  and  how  this  latter  becomes  a  function  of  just  the  shape  and  the  shape  controls 
because  of  the  nonholonomic  constraints.  We  show  that  the  configuration  and  shape  spaces  of 
the  (S'-snake  specify  a  principal  fiber  bundle  and  that  the  nonholonomic  constraints  specify  a 
connection  on  it. 

In  section  3  we  focus  on  2-module  G-snakes  (£  =  3)  evolving  on  3-dimensional  Lie  groups 
(n  =  3).  In  particular,  we  examine,  apart  from  SE( 2),  the  Heisenberg  group  H( 3),  the  Special 
Orthogonal  group  50(3)  and  the  Special  Linear  group  SL( 2).  We  derive  the  corresponding  Wei- 
Norman  representation,  the  system  kinematics,  the  connection  and  specify  the  nonholonomic 
singularities  in  each  case. 

In  section  4  we  discuss  possible  further  extensions  of  this  work.  For  reasons  having  to  do 
with  ease  of  exposition,  we  limit  ourselves  to  matrix  Lie  groups  in  this  paper.  Extensions  to 
arbitrary  Lie  groups  are  easy. 

2  Nonholonomic  Kinematic  Chains  on  Lie  Groups 

In  section  2.1  we  discuss  the  Wei-Norman  representation  of  curves  in  n-dimensional  Lie 
groups  and  in  section  2.2  we  derive  the  kinematics  of  the  (£  —  l)-module  kinematic  chain.  In 
section  2.3  we  examine  the  geometric  structure  of  the  chain  kinematics  when '  nonholonomic 
constraints  are  present,  using  the  theory  of  connections  on  principal  fiber  bundles. 

2.1  The  Wei— Norman  Representation  of  Lie  Groups 

Consider  a  left-invariant  dynamical  system  on  a  matrix  Lie  group  G  with  n-dimensional  Lie 
algebra  Q.  For  g  G  G,  the  left  translation  by  g  is  defined  as  the  map  Lg  :  G  — ►  G  :  h  t-»  gh,  for 
h  G  G.  If  e  is  the  identity  of  G,  then  TeLg  is  the  tangent  of  the  map  Lg  at  e.  Consider  a  curve 
g(.)  C  G.  Then,  there  exists  a  curve  £(.)  6  Q  such  that: 

9  =  TeLg  •  £  =  g  ■  £  .  (1) 
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Let  {Ai,  i  =  1, . . n)  be  a  basis  of  Q  and  let  [., .]  be  the  usual  Lie  bracket  on  Q  defined  by: 
[Ai,Aj\  =  AiAj  -  AjAi.  Then,  there  exist  constants  T*j,  called  structure  constants,  such  that: 

n 

[Ai,Aj]  =  ,  i,j  =  1, . . n.  (2) 

fc= i 


Let  Q*  be  the  dual  space  of  Q,  i.e.  the  space  of  linear  functions  from  Q  to  M.  Let  {A\,  i  = 
1, . . . ,  n)  be  the  basis  of  Q *  such  that 

A\{Aj)  =  6j,  i,j  =  1, . .  .,n,  (3) 

where  6j  is  the  Kronecker  symbol.  Then  the  curve  £(.)  C  Q  can  be  represented  as: 

n  n 

j=i  t=i 


Proposition  2.1.1  (Wei  k  Norman  [1964]) 

Let  5(0)  =  e,  the  identity  of  G  and  let  g(t)  be  the  solution  of  (1).  Then,  locally  around 
t  —  0,  g  is  of  the  form: 

g(t )  =  e^M1  e~/2(t)A 2  .  .  .  e-yn(t)An  ? 

where  the  coefficients  7 *  are  determined  by  differentiating  (5)  and  using  (1).  Then: 


f  71  N 
\7  nj 


(* A 


M(7i,...,7n) 


(6) 


The  matrix  M  is  analytic  in  7  and  depends  only  on  the  Lie  algebra  Q  and  its  structure  constants 
in  the  given  basis.  If  Q  is  solvable,  then  there  exists  a  basis  of  Q  and  an  ordering  of  this  basis, 
for  which  (5)  is  global.  Then  the  7;’s  can  be  found  by  quadratures. 


For  g  €  G,  £  €  G,  define  the  adjoint  action  of  G  on  Q  denoted  Adg  :  Q  — *•  Q  by: 

A dg{  d=  gtg'1  .  (7) 

From  (4)  we  have: 

3  3 

Ad ^  =  E  Abi(OAdgAi  .  (8) 

t=i  i=i 
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From  (5)  we  have: 


AdgAi  =  gAig-1  =  e71-41  •  •  ■e^AnAie-^A"  ■  ..e~^  . 


(9) 


Equation  (9)  can  be  made  more  explicit  by  the  Baker-Campbell-Hausdorff  formula  (Wei  & 
Norman  [1964]),  which  for  x,y  G  G  states  that: 


exye  x  =  y  +  [x,  y\  +  —[x,  [x,  y]]  +  dots 


(10) 


2.2  The  ( £  —  l)-module  Kinematic  Chain 

We  consider  a  dynamical  system  that  evolves  on  the  Cartesian  product  Q  —  G  x  •  •  •  x  G 

I  times 

Its  trajectory  is  a  curve  g(.)  =  (g^.),  ■ .  .,y  (.))  C  Q.  On  each  copy  of  G,  the  system  traces  a 
curve  <7.(.)  C  G ,  such  that 


9  —  TeLg  •  £  —  g  •  £_ ,  i  —  1, . . . ,  * 


(11) 


where  £.(.)  €  G,  i  =  1,  ...,£.  We  think  of  the  g's  as  the  nodes  of  a  kinematic  chain. 

X  1 

Let  the  instantaneous  shape  of  the  kinematic  chain  be  given  by  the  (£  —  l)~tuple 
i  f)  £  G  x  •••xG,  where 

(l— 1)  times 


g  • . .  _  g .  g  ■ , .  ? 

*,i+i  i  j+i 


(12) 


A  pair  of  adjacent  nodes  of  the  chain  constitutes  a  module.  The  g  ’s  can  be  regarded  as  the 

shapes  of  the  modules  of  an  (£  —  l)-module  kinematic  chain.  We  consider  the  corresponding 

curves  £  C  G,  specified  by  (1)  as  the  controls  of  the  system  (11).  We  refer  to  them  as  the 
i,i+l 

shape  controls: 


9 .  ■ , —  TeLg 


i,i+l 


i,t+l 


■  £  —9  •  £. . 

t,i+l  i,«+l  »,i+l 


*  =  !,...,£-  1 


(13) 


For  future  reference,  define  also: 


9.  .  =  9.  Y  =  9.  ...  ■■•9,  .  .i  *  <  j 

t  3  M+l  3-  1,3 


(14) 


and 


g.  .  =  TeL 


9.  . 


z i,j  -  9i'j 


i  <  j 


(15) 
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We  can  think  of  the  £ ’s  as  characterizing  the  global  motion  of  the  G-snake  system  with  respect 
to  some  global  coordinate  system,  while  the  £.  ,’s  capture  the  relative  motion  (or  shape  variation) 

*  ij 

of  nodes  i  and  j. 

From  (11),  (12)  and  (13)  we  get: 


f.  —  £.  .  .  +  Ad,  -i  , 

t  1,* 


-if  ,  i  =  2,... 


t-i 


(16) 


Applying  (16)  iteratively  we  can  express  any  as  a  function  of  £  and  of  the  shape  controls 

i  1 

£  ,...,£  as  follows: 

sl,2’ 


+  +'‘'+  ^1)2  +  •  (17) 

Using  (4)  we  get: 


J=1  J=1  8-1  8 

j=l  ’2*  j=i  1  8 


(18) 


2.3  Nonholonomic  Constraints  and  Connections  on  Principal  Fiber  Bundles 

In  this  section  we  consider  nonholonomic  constraints  acting  on  the  G-snake  and  we  show 
that  they  specify  a  connection  on  the  principal  fiber  bundle  associated  to  our  problem. 
Codimension  1  Constraint  Hypothesis:  Assume  that  the  evolution  of  system  (11)  on  each 
copy  of  G  is  constrained  to  lie  on  an  (n  —  l)-dimensional  subspace  h  of  the  Lie  algebra  Q,  where 
h  is  not  a  subalgebra  of  Q,  i.e. 

Z.eh,i  =  l, (19) 
Then,  for  some  AbK  €  Q *  (not  necessarily  an  element  of  the  basis  {A\,  i  =  1, . . .,  n})  we  have: 

h  =  Ker(^).  (20) 

The  constraints  (19)  can  be  expressed  as: 

4(£.)  =  0,  i  =  (21) 
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The  constraints  (21)  are  linear  in  the  components  of  £  and  those  of  the  shape  controls  £  . 

1  *>1 

This  can  be  made  explicit  by  defining  the  composite  velocity  vector  of  the  kinematic  chain: 


- d=  ( £x  ...e 

1  n  1  n 


ct-i,e  nt 


=  ( X(«,)  •  •  •  XU,)  X(fu)  •  •  •  X(ft_w) ) ' 


Proposition  2.3.1 

The  l  nonholonomic  constraints  (21)  can  be  written  in  matrix  form  as: 


(22) 


where  A  is  a  function  of  only  the  shape  of  the  system  and  is  a  block  lower  triangular  i  X  nl 
matrix  of  maximal  rank  of  the  form 


A  = 


/  *1,1 

0 

0  0  0 

0 

0  \ 

*1,2 

*2,2 

0  0  0 

0 

0 

0 

0 

0 

*1,» 

*2  ,i 

•••  *t-l,i  *»,» 

0 

0 

, 

... 

0 

\  *i,r 

*2,e 

•••  *»-i,  t  *f  ,e 

•••  *1-1,1 

*t,tj 

with  the  1  x  n  block  *p>q,  defined  for  p  <  q  as: 

*P,  q  =  (  ^It  (Ad^g-lg  )—  1  -^-1  ) 
\  Pi 


AK(Ad(g-ig  )-i^4„) 

p  ? 


(23) 


Proof 

From  (18)  and  (21)  : 


A^t)  =  '£A)(^)AbK(Aj)  =  0, 

i=i 

1= i  1=1  1 

n  n 

H - f  y;^(£i  n)A\{Ad{Q-lq  )-i^4j)  +  y Vi(£1Ml(Mr1».)-i^l)i 

j=l  ’  2  ’  1=1  1  1 

i  =  2, 


(24) 
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The  diagonal  blocks  *PiP  of  A  have  the  form  ^  ^(^4i)  •  •  •  AbK(An)  ^ ,  therefore  they  contain 

at  least  one  non-zero  constant  term.  Thus  A  has  always  maximal  rank. 

I 


Proposition  2.3.2 

Assume  £  >  n.  Partition  E  as  (Ei  E2),  with  E2  an  ^-dimensional  vector  containing  the 
components  of  £  (and  possibly  some  components  of  shape  controls),  while  Ei  is  an  (n  —  1)1- 
dimensional  vector  containing  only  components  of  shape  controls.  Let  the  corresponding  partition 
of  A  be  (Ai  A2),  with  A\  a  (n  -  1)^  x  l  matrix  and  A2  a  locally  invertible  l  x  £  matrix.  Then 
from  (22)  : 

“2—  -d-2  (si2>  •  •  • 

Proof 

Follows  from  the  smooth  dependence  of  A  on  the  shape  variables  and  the  maximal  rank  property 
of  Proposition  2.3.1. 


The  physical  significance  of  this  result  is  that,  if  the  global  motion  of  the  {£  —  l)-module 
system  is  characterized  by  the  global  motion  of  its  first  module  (i.e.  by  f  ),  then  variations  of  the 
shape  controls  (at  least  those  which  are  elements  of  Ei)  induce  a  global  motion  of  the  system. 
Configurations  where  A2  is  singular  will  be  called  nonholonomic  singularities. 

Definition  2.3.3  (Principal  Fiber  Bundle  (Nomizu  [1956])  ) 

Let  5  be  a  differentiable  manifold  and  G  a  Lie  group.  A  differentiable  manifold  Q  is  called 
a  (differentiable)  principal  fiber  bundle  if  the  following  conditions  are  satisfied: 

1)  G  acts  on  Q  to  the  left,  freely  and  differentiably: 

$  :  GxQ  -*Q  :  {g,q)  i->  g  •  qd=  •  q  . 

2)  S  is  the  quotient  space  of  Q  by  the  equivalence  relation  induced  by  G,  i.e.  S  =  Q/G  and  the 
canonical  projection  7r  :  Q  — ►  S  is  differentiable. 

3)  Q  is  locally  trivial,  i.e.  every  point  s  €  S  has  a  neighborhood  U  such  that  7T— 1(f7)  C  Q  is 
isomorphic  with  U  x  G,  in  the  sense  that  q  €  7T-1(f7)  i->  (^(g),  <j>(q))  €  U  x  G  is  a  diffeomorphism 
such  that  4> :  7t_1(17)  — +  G  satisfies  <f)(g  •  q)  =  g<p{q)^g  G  G. 


8 


For  s  £  S,  the  fiber  over  s  is  a  closed  submanifold  of  Q  which  is  differentiably  isomorphic  with 
G.  For  any  q  £  Q,  the  fiber  through  q  is  the  fiber  over  s  =  x(<?).  When  Q  =  S  X  G,  then  Q  is  said 
to  be  a  trivial  principal  fiber  bundle. 

I 

Consider  now  the  manifolds  Q  and  S  defined  in  section  2.2  and  the  canonical  projection 
x  :  Q  — s ►  S  defined  by  equation  (12),  i.e. 


(26) 


We  can  show  that  the  quadruple  (Q,  S,n,G),  together  with  the  action  $  of  G  on  Q  defined  by 

$  ;  GxQ  -+Q  :  ( g,q)=  {g,  (g^ . . .  ,g{))  »-*■  g-q=  (gg^. .  .,gg£)  ,  (27) 

meets  the  requirements  of  Definition  2.3.3,  and  therefore  ( Q,S,n,G )  is  a  trivial  principal  fiber 
bundle.  Indeed,  x  is  differentiable  and  its  differential  is 


7I’*«  '  TqQ  :  (fl'1£1) 


^1,2  ^1,2’ 


where  the  £.  i  .  are  given  by  (16)  : 


(28) 


,  ,  =  f.  ~  Ma-l 


i— l,i 


(s-1  g.y 
2  —  1 


i  =  2, . .  .,£  . 


(29) 


Definition  2.3.4  (Connection  on  a  Principal  Fiber  Bundle  (Nomizu  [1956])  ) 

Let  (Q,S,  x,G)  be  a  principal  fiber  bundle.  A  connection  on  the  principal  fiber  bundle  is 
a  choice  of  a  tangent  subspace  Hq  C  TqQ  at  each  point  q  £  Q  (horizontal  subspace)  such  that, 

if  Vq  d=  {u  £  TqQ |x*4(u)  =  0}  is  the  subspace  of  TqQ  tangent  to  the  fiber  through  q  (vertical 
subspace),  we  have: 

1)  TqQ  =  Hq  ©  Vq. 

2)  For  every  g  £  G  and  q  £  Q,  Tq$g  ■  Hq  =  Hg.q. 

3)  Hq  depends  differentiably  on  q. 
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Proposition  2.3.5 

Away  from  nonholonomic  singularities  and  when  £  =  n,  the  nonholonomic  constraints  (21) 
specify  a  connection  on  the  principal  fiber  bundle  (Q,  S,  tt,  G),  with  the  horizontal  subspace 
defined  as  follows: 


Hq  =  {veTqQ  |  t>  =  and  f.  6  h } 


{veTqQ\v  =  (g^x,..^9e^)  and  S2  =  -A2\Tr(q))A1(n(q))E1  }  , 


(30) 


where  =  (^-2  •  •  •  ^-2  £2-3  •  •  •  ^)T  and  S2  =  (f1  •  •  •  ^)T  . 

1  n  1  n  1  n 

Proof 

Due  to  the  left-invariance  of  our  system,  TqQ  =  {(<7  £  , . .  -,gg£g)  I  £.  £  £}.  The  vertical  subspace 
is  (from  (27)  —  (29)) 


v,  =  {V  <S  I  TT,^  =  0} 

—  {(9^1  ■  •  1  (^1)2^12’  •  •  •’  y 

(31) 

—  {(^^j >  •  •  •  ’  I  ^1>2  —  —  ^e-itt  ~  ^ 

=  {(51C1?  •  •  '’9(te)  I  =  Ad^g-.ig^-i£,x,  i  =  2,  ...,£}  . 

Physically,  the  vertical  subspace  contains  all  infinitesimal  motions  of  the  kinematic  chain  that 
do  not  alter  its  shape. 

To  show  property  (1)  of  Definition  2.3.4,  we  first  prove  that  Hq  fl  Vq  =  {0}  and  then  that 
dim(T,Q)  =  dim  (Hq)  +  dim(Vg). 

To  show  Hq  nFg  =  {0},  assume  that  there  exists  a  non-trivial  v  =  q-£^Hqn  Vq.  By  the  definition 
of  Vq,  the  corresponding  shape  controls  are  zero.  Thus  Si  =  0  and,  by  the  definition  of  Hq,  also 
~2  =  0.  But  then  £  =  0  and  from  (31)  also  £.  =  0,  i  =  2, . .  Thus  £_0.  Thus  Hq  f)Vq  =  {0}. 
Now  observe  that,  away  from  the  nonholonomic  singularities  dim(ifg)  =  n£  —  £  .  Further, 
dim(Vg)  =  n  .  So,  when  £  =  n,  dim(ffg  @  V,)  =  ( n£  —  £)  +  n  —  (n2  —  n)  +  n  =  n2  =  dim(TgQ).  It 
follows  that  Hq  ®  V,  =  TqQ. 

To  show  property  (2)  of  Definition  2.3.4,  consider  Tq$g  •  Hq  =  g  •  Hq  =  g  ■  {(g^^  ■  ■  - ,  ge£e)  \  € 

h}  d=  {(gg^x,-  ■  -,gge^)  I  t  e  h}  and  Hg.q  =  {v  €  Tg.qQ  \  v  =  (g  •  q)  •  (^,...,^)  and  £.  € 
h }  =  {{gg^-.^gg^)  If.  G  h}  .  Then,  obviously,  Tq$g-Hq  =  Hg.q. 
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Property  (3)  of  Definition  2.3.4  is  immediate  from  the  smooth  dependence  of  A  on  the  shape  and 
from  the  left-invariance  of  our  system. 


3  Three-dimensional  Nonholonomic  Lie  Groups 

Here  we  specialize  the  results  of  the  previous  section  to  kinematic  chains  on  Lie  Groups 
with  3-dimensional  real  non-abelian  Lie  algebras  (n  =  3).  In  section  3.1  we  consider  the  Special 
Euclidean  group  SE( 2),  in  section  3.2  the  Heisenberg  group  H( 3),  in  section  3.3  the  Special 
Orthogonal  group  S0( 3)  and  in  section  3.4  the  Special  Linear  group  SL( 2). 

We  study  2-module  (£  =  n  =  3)  kinematic  chains  on  each  of  these  groups  by  deriving 
their  Wei-Norman  representation  and  by  defining  the  partial  connection  on  the  corresponding 
principal  fiber  bundle. 

Let  G  be  one  of  the  above  four  matrix  Lie  groups  and  Q  be  the  corresponding  Lie  algebra. 
Consider  the  system  (2.1)  on  G  : 

9  —  TeLg  •  £  =  g  •  £  ,  (1) 

with  g  €  G  and  £  €  Q.  From  Proposition  2.1.1,  any  g  6  G  has  a  local  Wei-Norman  representation 
of  the  form 

g(t)  =  e7i(<Mie72(<M2e73(iM3  _  ( 2 ) 


Consider  now  the  2-module  kinematic  chain  on  G.  From  the  system  kinematics  (equations 
(2.11)  —  (2.18))  we  have: 


^2  ^1^1,2  ’ 


0=00  =  a  g  a 

a3  a2  2,3  alal,2  2,3 


(3) 


g  —  g  g 

i,3  i,2J2,3 

From  (2.16)  we  get  for  the  corresponding  velocities: 


C  £.  o  +  Ad{g-ig  \  i£  , 


1,2 


'l  *2' 


^3  ^2,3  +  AdG~'9z)  ’2  ’2,3 


"*£.  “£.<.+  Ad{g-'gA~'t i  2  +  • 


(4) 


2  V  ’1.2  "3' 

Assume  that  the  evolution  of  system  (1)  on  each  copy  of  G  is  constrained  to  lie  on  a  2-dimensional 
subspace  h  of  the  Lie  algebra  Q ,  where  h  is  not  a  subalgebra  of  Q.  Define 

=i  =  (e2  e2  e2  e3  e3  e3)T  and  ~2  =  &  e  ey . 

1  i  O  1  <i  o  1^0 
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Proposition  2.3.1  holds  with  E  =  (Ej  E2).  From  Proposition  2.3.2  we  conclude  that  the  global 

velocity  of  the  2-module  kinematic  chain,  as  it  is  characterized  by  £  ,  can  be  expressed  as  a 

function  of  only  the  shape  variables  g  ,g  and  shape  controls  £  ,  £  of  the  assembly: 

1)2  2,3  1 ,2  2,3 

“2  —  A2  (9l  2i92  3)-^1(9i  2^S2  3)—i  •  (5) 

From  Proposition  2.3.5,  equation  (5)  defines  (away  from  the  singularities  of  A2)  a  connection  on 
the  trivial  principal  bundle  (5  X  G,S,n,G)  with  S  =  G  x  G  and  with  horizontal  subspace: 


Hq  =  {v  6  TqQ  |  «  =  (g^g^g^)  and  £.  €  h} 

=  {veTqQ  \v  =  (g^g^g^)  andS2  =  -A2l{g^g^)  } 


(6) 


Subsequently,  we  will  derive  explicitly  the  Wei-Norman  representation  for  each  of  the  Lie 
groups  mentioned  earlier  and  we  will  define  the  connection  (6)  for  specific  2-dimensional  sub¬ 
spaces  h  of  Q.  In  (Vershik  &  Gershkovich  [1994])  the  authors  present  a  result  showing  that, 
for  each  of  H (3),  50(3)  and  SE( 2),  all  2-dimensional  subspaces  h  of  Q,  which  are  not  subal¬ 
gebras,  are  isomorphic  and  can  be  represented  by  h  —  sp{.4i,.42}  in  the  basis  of  Q  specified 
in  the  following  sections.  For  SL( 2),  there  are  2  classes  of  such  isomorphic  subspaces  that  can 
be  represented,  respectively,  by  h  =  sp{^4i,^2)  and  by  h  =  sp{yl3,^4i  +  ^2}-  Therefore,  only 
nonholonomic  constraints  corresponding  to  these  subspaces  of  Q  will  be  considered  here. 

Our  main  purpose  in  this  section  is  to  set  the  stage  for  a  deeper  understanding  of  this 
novel  class  of  kinematic  chains,  by  cataloguing  the  low-dimensional  possibilities.  One  case, 
corresponding  to  SE( 2)  has  already  found  a  concrete  mechanical  realization  (Krishnaprasad 
&  Tsakiris  [1994]).  Others  might  follow,  for  instance,  there  are  possible  connections  between 
50(3)-snakes  and  the  kinematics  of  long  chain  molecules  (Karplus  &  McCammon  [1986]). 

3.1  Nonholonomic  Kinematic  Chains  on  the  Special  Euclidean  Group  SE(2) 

Let  G  —  SE( 2)  be  the  Special  Euclidean  group  of  rigid  motions  on  the  plane  and  Q  =  se{ 2) 
be  the  corresponding  algebra  with  the  following  basis: 

/  0  -1  0\  /O  0  1\  /O  0  0\ 

^1=1  0  0  ,  A2  =  0  0  0  ,  ^3  =001.  (7) 

\o  0  0/  \0  0  0 )  \0  0  0/ 


Then: 


[Ai  ,  A2]  =  ^43  ,  [Ai  ,  ^43]  —  —  *42  )  [A2  ,  ^3]  —  0  . 


(8) 
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The  algebra  Q  is  solvable  and,  from  Proposition  2.1.1,  any  g  6  G  has  a  global  Wei-Norman 
representation  of  the  form  (2),  with 


71 

72 

73 


Equation  (9)  is  solvable  by  quadratures: 


(9) 


7i0)  =  7i(0)  + 


5 


t  t 

72 (t)  =  72(0)  cos  ^  J  ^(o-)d(T^  +  73(0)  sin  ^  J  a)da 


+ 


t  t  t  t 

J  £20)cos^y  (i(<T)dajdT  +  J  £3(  r)  sin  (^J^)d(^jdr,  (10) 


t  t 

73 (0  =  -72(0)  sin  (^J  ^(cr)do^  +73(0)  cos  (^J  ^(c r)da 


t  t  t  t 

-  J  ^)sin(/ ^(o)dajdr  +  J  ^(t)cos^J  ^(a)dajdr 


From  (2.9),  (9)  and  (10)  we  can  compute  Adg-iAi  : 


Adg-iAi  =  A\  —  73^2  +  72*43  > 

Adg-iA2  =  cos 71^2  -  sin  7i*43  ,  (11) 

Adg-xAz  =  sin7i*42  +  cos 71^3  . 


Consider  now  the  2-module  kinematic  chain  on  SE( 2),  (£  -  3).  A  concrete  mechanical  realization 
of  such  a  system  is  the  Variable  Geometry  Truss  assembly  mentioned  in  section  1,  which  appears 
in  fig.  3.1.  The  system  kinematics  of  equations  (3),  (4)  apply. 

From  (8)  we  can  see  that  there  are  two  possible  2-dimensional  subspaces  h  of  Q  that  can  generate 
the  whole  algebra  under  Lie  bracketing: 

h3  -  sp {A\ , .A  2 }  =  Ker(^3)  and  h2  =  sp{*4i,^43}  =  Ker^)  •  (12) 
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Fig.  3.1 


Subsequently  we  will  consider  only  h2  (which  is  exactly  the  case  of  the  system  in  fig.  3.1).  The 
nonholonomic  constraints  £.  €  h2  can,  then,  be  expressed  as: 


4(0  =  0,  *=1,2,3. 


Equation  (5)  holds  with  Ex  and  E2  defined  as  above  and  with: 


Ax  = 


0  0  0' 
0  0  0 


0  0  0 

0  1  0 

4 (Ad(3-i3  )-i Ax)  4 (Ad(5-is  )-i A2)  A\(Ad(g-ig  yiAz)  0  10 

1  2  3  2  3  2  3 

0  0  0  0  0  o' 

0  1  0  0  0  0 

—J3'3  cos7j’3  sin7^’3  0  10 


and 


A.2  — 


\ 


0  1  0 

Al(Ad(g-igj-iAi)  4(Ad(ff-ia2)-iA2)  4(Ad(5-i52)-iA3) 

yA\{Adlg-lg j-iAx)  4 (Ad(_i)-i.A2)  4(Ad(  -i  )-iA3)  ^ 

'  1  J  1  O  1  d 


1 


0 

sin  7  { 


3  3  cos\, 

73  ’  c°s7i"  sm  7{ 


3  _ :,1.3 


(13) 


The  nonholonomic  singularities  of  the  system  are  the  configurations  where: 

det(A2)  =  4 {Ad(g-ig  )-1Ai)AL2{Ad^g--Lg  )-iA3)  -  A\{Ad^g-ig  yiA\)A\{Ad^g-ig  yi.Az) 

13  12  12  1  J 


=  ~73 ’3  sin  74  +  73’z  sin  7* 


1,2 


1,2 

1/  ’  i 


1,3 
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See  (Krishnaprasad  &  Tsakiris  [1994])  for  the  mechanical  interpretation  of  the  nonholonomic 
singularities.  Away  from  those,  equation  (6)  specifies  the  connection  corresponding  to  this  sys- 
tem.  Under  periodic  shape  controls  £12>€2  3)  the  system  translates  when  ’  =0  and  rotates 

otherwise.  This  in  turn  is  associated  with  the  concept  of  geometric  phase  of  the  system,  which 
is  discussed  in  greater  detail  in  (Krishnaprasad  &  Tsakiris  [1994]).  In  the  case  of  the  translating 
system,  the  corresponding  geometric  phase  is  plotted  in  fig.  3.2  for  one  period  of  the  shape 
controls  and  for  varying  amplitudes. 


Geometric  Phase:  Translation 


1 


Fig.  3.2 


3.2  Nonholonomic  Kinematic  Chains  on  the  Heisenberg  Group  H(3) 

Let  G  =  H(3)  be  the  Heisenberg  group  of  real  3x3  upper  triangular  matrices  of  the  form 

1  a  /3 

0  1  7  ]  and  let  Q  =  h( 3)  be  its  algebra  with  the  following  basis: 

0  0  1 


0  1  0\  /  0  0  0 

Ai=[0  0  0  ,  Ai  =  0  0  1 

000/  \0  0  0 


0  0  1 
A3  =  I  0  0  0 
0  0  0 


(14) 


Then: 

[Ai  ,  A2]  =  A3  ,  [A\  ,  A3\  =  0  ,  [ A2  ,  A3]  =  0  .  (15) 

The  algebra  Q  is  nilpotent  (thus  solvable)  and,  from  Proposition  2.1.1,  any  g  €  G  has  a  global 
Wei-Norman  representation  of  the  form  (2)  with 


71  \  /I  0  0 

72  =  0  10 

,73/  \-72  0  1 


(16) 
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Equation  (16)  is  solvable  by  quadratures: 


t 

7i(0  =  7i(0)  +  J^r)dT  , 
o 

t 

72 (<)  =  72(0)  +  J  $2(r)dr  , 

0 

t  t 

7s(0  =  73(0)  -  J  72 (r)^(r)dr  +  J  f3(r)dr 

0  0 

t  t  T  t 

=  73(0) -72(0)  J  Zx{?)dcT-  j  ^)(/ ^(cr)dcrjdr+  J  ^(r)dr  . 

0  00  0 

From  (2.9),  (16)  and  (17)  we  can  compute  : 

Ad3-i.4i  =  Al  +  72^3  > 

Adg-iAi  —  A2  —  71^3  5 
Adg-iAz  =  ^3  . 


(17) 


(18) 


Consider  now  the  2-module  kinematic  chain  on  G,  (£  =  3).  The  system  kinematics  of  equations 
(3),  (4)  apply. 

From  (15)  we  can  see  that  there  is  only  one  possible  2-dimensional  subspace  h  of  Q  that  can 
generate  the  whole  algebra  under  Lie  bracketing: 


h  =  sp{>li,^42}  =  Ker^g)  . 

The  nonholonomic  constraints  can,  then,  be  expressed  as: 

4(^  =  0,  *  =  1,2,3. 

Equation  (5)  holds  with  £1  and  E2  defined  as  above  and  with: 


(19) 


(20) 


Ai 


/  0  0  0  0  0  0\ 

I  0  0  1  0  0  0  ] 

yAb3(Ad(g-lg)-iAi)  Aj  {Ad(g-lgyiA2)  A3(Ad(g-ig  )-i-43)  0  0  1 J 

2  2  o  2  o 


0  0  0  0  0\ 

0  10  0  0 
—7 1’3  10  0  1/ 
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and 


A2  = 


0  0  1  \ 
)-^i)  ^(Ad^g-ig J-1A2)  Ab3(Ad,{g-lgj-iA3) 


T  2 


"1  v 


1  2 


yA3(Ad(3_i3  )— i^4.i)  A3(Ad(3_i3  )-iA2)  -43(Ad(3_l3  yiA3)  / 

1  3  1  3  1  3 


0 


0 


1 


I  1,2  1,2  1 

7l3  -71,3 

.V  — 7i  1. 

The  nonholonomic  singularities  of  the  system  are  the  configurations  where: 

det(A2)  =  A^Ad^-i^  )-iAi)A3(Ad(3_i3  )-iA2)  -  A^Ad^-i^  ^-iAi)Ah3(Ad^g-lg  )-iA2) 
1  2  1  3  1  3  1  2 


1,2  1,3  ,  1,2  1,3  1,2  2,3  ,  1,2  2,3  n 

=  -72  7i  +7i  72  =-72  7i  +7i  72  =°- 


3.3  Nonholonomic  Kinematic  Chains  on  the  Special  Orthogonal  Group  SO(3) 

Let  G  =  £0(3)  be  the  Special  Orthogonal  group  of  real  orthogonal  3x3  matrices  with 
determinant  equal  to  one  and  let  Q  =  so(3)  be  the  algebra  of  3  x  3  real  skew-symmetric  matrices. 
Consider  the  following  basis  for  Q  : 


/  0  0  0  \ 

Ai  =  0  0  -1  ,  A2 

\0  1  0  ) 


(  0  0  1\  /O  -1  0\ 

0  o  o  ,  a3  =  1  0  0 

\-l  00/  \0  0  0 / 


Then: 


[Ai  ,  A2]  =  A3  ,  [Ai  ,  A3]  =  — A2  ,  [A2  ,  A3]  =  Ai  . 


(21) 


(22) 


The  algebra  Q  is  simple,  thus,  the  Wei-Norman  representation  (2)  is  only  local  (defined  when 
cos  72  i=-  0)  with  coefficients: 


(7i  \  /  sec 72  cos 73  -sec72sin73  0 \  /  ^  \ 

72  =  sin73  COS7 3  0  )  .  (23) 

73  /  \  -  tan  72  cos  73  tan  72  sin  73  1  /  \  £  / 

From  (2.9)  and  (23)  we  can  compute  Adg-iAi  : 

Adg-iAi  =  cos  72  cos73Ai  -  cos  72  sin73A2  +  sin72A3  , 

Adg-\A2  =  (sin7i  sin  72  cos73  +  cos  71  sin73)Ai 

+  (—  sin7i  sin  72  sin  73  +  cos  71  cos73)A2  —  sin  71  cos72A3  ,  (24) 

Adg-iA3  =  (-  cos 71  sin  72  cos 73  +  sin  71  sin 73)Ai 

+  (cos  7x  sin  72  sin  73  +  sin  71  cos  7 3)A2  +  cos  71  cos  72  A3  . 
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Consider  now  the  2-module  kinematic  chain  on  G,  (£  =  3).  The  system  kinematics  of  equations 
(3),  (4)  apply. 

From  (22)  we  can  see  that  there  are  three  possible  2-dimensional  subspace  h  of  Q  that  can 
generate  the  whole  algebra  under  Lie  bracketing: 


h3  =  sp{Ai,A2}  =  Ker^g),  h2  -  sp{^!,^3}  =  Ker(^), 
and  hi  =  sp{A2,-43}  =  Ker(_43), 


(25) 


We  consider  only  h3  C  Q ■  The  nonholonomic  constraints  f .  £  h3  can,  then,  be  expressed  as: 


4(y  =  0,  *  =  1,2,3. 


(26) 


Equation  (5)  holds  with  and  E2  defined  as  above  and  with: 


At  = 


0  0  0  0  0  O' 

0  0  1  0  0  0 

A^Ad^g-lg  )-l-4i)  A2(Ad^g-lg  )  ~  1  ^  ^  ^  3  (  A  g  -  1  g  ^  ~  1  ^  3  )  0  0  l 

2  3  2  3  2  3 


0 

0 


,sm  72 


2,3 


0 

0 

2,3, 


0  0  0' 
0  0  0 


sm7{’'  cos 72 ’3  cos7^’3  cos72’3  0  0  1 


and 


A2  — 


(  0  0  1 
A3{Ad(g-ig  j-iAi)  A2(Ad(g-igj-iA2)  Ab3(Ad{g-lg  j-iA3) 


T  “r 


T  2 


T  2 


^A3(Ad(s-i  g  )-iAi)  A3(Ad(g-ig  )-iA2)  A3(Ad(g~i g  s-iA3)  J 

1  O  1  1  O 


0 


0 


=  |  sin7,1,2  —  sin7j’2  COS72’2  cos7j’2  cos  79 ’2 
sin  72  ’3  —  sin  7j  ’3  cos  72  ’3  cos  7X  ’3  cos  72  ’3 


The  nonholonomic  singularities  of  the  system  are  the  configurations  where: 

det(A.2)  =  A3(Ad(g-ig  ^-iAi)A3(Ad^_lff  ^-1^2)  -  A3(Ad^g- lg  ^-iAi)A3(Ad^g-ig  )-i A2) 
1  2  1  3  1  3  1  2 


-  sin  72  ’2  sin  7^  ’ °  cos  72  ’ J  +  sin  72  ’ J  sin  7^  cos  72  =  0  . 


.,1,3  , 


1,3 


1-3, 


,1,2 


1,2 

V  ’  — 
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3.4  Nonholonomic  Kinematic  Chains  on  the  Special  Linear  Group  SL(2) 

Let  G  =  S L(2)  be  the  Special  Linear  group  of  real  2x2  matrices  with  determinant  one  and 
let  Q  =  5/(2)  be  the  algebra  of  real  2x2  matrices  of  trace  zero.  Consider  the  following  basis  for 


Q: 


(27) 


Then: 

[A\  ,  A2]  =  2A3  ,  [-4i  ,  .43]  =  —  .4i  ,  [A2  ?  -43]  =  -42  •  (28) 

The  Wei-Norman  representation  (2)  is  only  local  (defined  when  e73  7^  0)  with  coefficients: 


/  7i  \  (  ^  0  0\  nx\ 

72  =  72e73  e"73  0  U2  (29) 

Via/  \-272e73  0  1/  \HJ 

(See  however  comments  in  (Wei  &  Norman  [1964])  and  their  Theorem  3.  A  global  representation 
of  SL{ 2)  can  be  obtained  using  {-4i,  A\  —  ^2,^43}  as  a  basis). 

From  (2.9)  and  (29)  we  can  compute  Adg-\Ai  : 

Adg-iAi  =  e"7Mi  -  ~i2e'zA2  +  272A3  , 

Adg-iA2  =  ~7ie_73-4i  +  (7172  +  l)2e73-42  -  27^7172  +  l)-43  ,  (30) 

Adg-iAz  =  7iC_73Ai  -  72(7172  +  l)e73-42  +  (27x72  +  l)-43  . 

Consider  now  the  2-module  kinematic  chain  on  G,  (£  =  3).  The  system  kinematics  of  equations 
(3),  (4)  apply. 

From  (28)  we  can  see  that  there  are  two  possible  2-dimensional  subspaces  h  of  Q  that  can  generate 
the  whole  algebra  under  Lie  bracketing: 

h3  =  sp{Ax,-42}  =  Ker(4)  and  h\,2  =  sp{*43,-4i  +  -42)  =  Ker(4  -  4)  ■  (31) 

We  first  consider  h3  C  Q-  The  nonholonomic  constraints  f .  £  h3  can,  then,  be  expressed  as: 


4(^  =  0,  <  =  1,2,3. 

Equation  (5)  holds  with  Si  and  S2  defined  as  above  and  with: 

/  0  0  0  0  0  0\ 

A  _  j  0  0  1  0  0  0  I 

l  A\{Ad(g-ig  )-iAx)  4 {Ad^g-'g  )-1^)  4(^(s-15  )— ^3)  0  0  1 J 

'  2  3  2  3  2  3 


/  0  0  0  0  0  o\ 

j  0  0  1  000 

V2722’3  -27x2’3(7i’3722’3  +  1)  27i2’3722'3  +  1  0  0  1  J 


(32) 
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and 


0  0  1 
•Aa{Ad(g-igj-iAi)  Al(Ad(g-igj-iA2)  A^Ad^  g^-iA3) 

•AiAd^-!  )-iAi)  A$(Ad(g-i  yiAz)  Ab3(Ad^g-lg  )-i-43) 


4  3  -271  yhl  3 + 1}  2i'li'l  +  1 

272’  -27i’(7i’72’  +  1)  27l1’3721-3  +  1 


The  nonholonomic  singularities  of  the  system  are  the  configurations  where: 


det(^2)  =  Az{Ad(g-lg  )-iA\)A\(Ad(g-ig  )-iA2)  -  A^Ad^-t g  )-iA\)A\(Ad(g-1  g  )-'A2) 

1  2  1  3  13  1  * 


=  -4721’27i1,3(7i1'3721’3  +  1)  +  4721,37i1’2(7i1’2721’2  +  1)  =  0 


Now  we  consider  the  subspace  h\t2  C  Q  in  (31).  The  nonholonomic  constraints  can,  then,  be 
expressed  as: 


(4-4)(0  =  °,  *  =  1,2,3. 


Equation  (5)  holds  with  Si  and  E2  defined  as  above  and  with: 


A\  = 


0  0 

1  -1 

(A  -  A){A4(,-,.  ,-^i)  (4  -A\){Ad(t„  r,Ai) 

2  3  2  3 


0  0  0  0 

0  0  0  0 

{A\  -  A\)  (Ad(a_i g  )-'Az)  1-10 

2  3 


o  o 

-3’3  +  (722’3)2e^3  -‘(7112’3722’3  +  1)2<^ 


0  0  0  0 

0  0  0  0 

7i’V'1'3'3  +  722’3(7l2'3722'3  +  l)e^2'3  1  -1  0 
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and 


A2 


1  -1 

(A  -  A)  (AiU-'hr'A')  M  -  A) 

1  (A  -  A) (M(r^r'A)  (A  -  A)  (^(»-..  )-*■ A) 

0  \ 

(A\  -  A2)  (Ad{g-lg^-iA3) 

(■^1  -  A\)  (Ad(g-lg^-iA3)  j 


e  V  -f  (72’2)2e'73’ 

e-7^’3  -f  (72,3)2e73’ 


-1 


-(Ti1’2)2^'2  -  {ll'WS  +  l)2e73 

-(7i1>3)2e-^3  -  (7i1’3721’3  +  l)2^1 


7i1’2e-73'2+72’2(7i1’2721’2  +  l)^’2 
7j’3e-73'3  +  72'3(7i’372’3  +  l)^7® 3 


The  nonholonomic  singularities  of  the  system  are  the  configurations  where: 
det(A2)  =  (^4i  -  Ab2)(Ad(g-ig  )-i(Ai  +-42))(^4i  -  A\)  (Ad^g-ig  )-i^3) 

1  2  1  3 

-  (*4i  -  A2){Ad(g-ig  +  A2))(Ai  -  A2)(Ad(g-igj-iA3) 

1  3  1  2 


=  (<r’“  +  (ASM*  -  -  (tiAV + i)v:’) 


4  Conclusions 

In  this  paper  we  introduce  the  concept  of  G-snakes,  which  is  a  class  of  kinematic  chains 
with  nonholonomic  constraints  evolving  on  a  Lie  group.  Shape  variations  of  the  system  modules 
induce  a  snake-like  global  motion  of  the  system.  We  provide  the  framework  upon  which  motion 
planning  strategies  based  on  periodic  shape  variations  can  be  developed.  A  concrete  mechanical 
realization  is  associated  with  G  =  SE( 2).  We  offer  a  catalogue  of  low-dimensional  possibilities. 
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Further  extensions  of  this  work  include  the  study  of  geometric  phase  for  each  of  the  groups 

that  we  discuss,  as  well  as  the  study  of  optimal  control  problems  related  to  the  optimal  choice 

of  shape  variations  that  will  achieve  motion  between  two  desired  configurations. 
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